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ABSTRACT 

Massive protostars have associated bipolar outflows that can produce strong shocks when they interact with the surrounding medium. 
At these shocks particles can be accelerated up to relativistic energies. Relativistic electrons and protons can then produce gamma-ray 
emission, as some theoretical models predict. To identify young galactic objects that may emit gamma rays, we crossed the Fermi 
First Year Catalog with some catalogs of known massive young stellar objects (MYSOs), early type stars, and OB associations, and 
we implemented Monte Carlo simulations to find the probability of chance coincidences. We obtained a list of massive MYSOs that 
are spatially coincident with Fermi sources. Our results indicate that ~ 70% of these candidates should be gamma-ray sources with a 
confidence of ~ 5<x. We studied the coincidences one by one to check the viability of these young sources as potential counterparts to 
Fermi sources and made a short list of best targets for new detailed multifrequency observations. The results for other type of young 
galactic objects are not conclusive. 
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1. Introduction 

Recently, massive young stellar objects (MYSOs) have been 
suggested as gamma-ray sources (|Araudo et aL|2007| |Romero| 
2008, |Bosch-Ramon et aL]|2010^ Massive stars are formed 
in dense cores of cold clouds. The processes that take place 
during the formation of the star are mostly unknown. It is 
clear, however, that the formation of massive stars involves 



outflows (Gar ay & Lizano|[T999| |Reipurth & BaUypO Ol). The 
accumulation of material around the core of the cloud would 
generate a massive protostar that starts to accrete material from 
the environment. The accretion is expected to have angular 
momentum that leads to the formation of an accretion disk. 
The rotation would twist the strong magnetic fields present 
in the progenitor cloud around the disk, where a magnetic 
tower can be formed, giving rise to collimated outflows or 
jets, as simulations predict (Banerjee & Pudritz 2006 2007). 
The observational evidence of outflows comes from methanol 
masers and from direct detection of thermal radio jets. These 
jets propagate along distances in a fraction of a parsec (Marti 
|et al.|[T993] ). At the jet termination region, interaction with the 
external medium creates two shocks: a bow shock moving in the 
interstellar medium (ISM) and a reverse shock in the jet. These 
shocks can accelerate particles that, in turn, can produce gamma 
rays trough inverse Compton (IC) scattering of infrared (IR) 
photons, relativistic Bremsstrahlung or inelastic proton-proton 
collisions, if protons are accelerated as well. In some cases 
non thermal radio lobes and jets have been observed, indicating 
the presence of relativistic electrons that produce synchrotron 
radiation ( Garay et al.|2003 Carrasco-Gonzalez et al.| 2010). 



Other possible scenarios have been suggested for the 
gamma-ray production involving young stars, such as the case 



of the massive stars with strong winds. In this scenario, gamma- 
rays could be produced in the interaction between the supersonic 
winds and the ISM. The terminal shock can accelerate particles 
and ions up to high energies, which might interact with the ambi- 
ent matter producing gamma-rays. Whereas the luminosity pro- 
duced by a single massive star wind should be low, collective 
effects might be important ( Torres et al.|2004 1. 
Gamma-rays can also be produced in the wind interaction region 
of a WR+OB binary system (Benaglia & Romero|2003| l. In this 
case, the acceleration region is between the two components of 
the binary system and is exposed to strong photon fields where 
IC cooling of the electrons can generate a significant amount of 
high-energy (HE) non thermal emission. A source of this class, 
77 - Carinae, has been recently detected at £ > 100 MeV by 
AGILE ( |Tavani et al.|2009) . 

Of-type stars have strong winds with velocities higher 
than the escape velocity, which implies a strong mass loss rate 
(10~ 6 - 10~ 5 M Q yr ). The action of this wind in the interstellar 
medium can create hot gas bubbles with expanding boundaries 
of swept-up material, which might produce gamma rays in a 
similar way to the case of WR stars. The case of gamma-ray 
emission in Of-type stars has been discussed in the past, e.g. by 
Voel k & Forman| (|l982). The predicted luminosity, however, is 
still below the current sensitivity of gamma-ray instruments. 

Finally, OB associations are tracers of a number of galactic 
objects that can produce gamma rays, such as neutron stars, 
massive stars with strong winds, young stellar objects (YSOs), 
etc. They are also thought to be places where acceleration of a 
significant fraction of galactic cosmic rays (CRs) might occur 
(e.g. |Binns et al.|2008) . 
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The aim of this work is to find evidence supporting the pres- 
ence of HE emission coming from massive YSOs and other 
young galactic sources. To attain this goal we study the spa- 
tial coincidence between gamma-ray sources detected by Fermi 
and samples of young objects, such as YSOs, WR stars, Of- 
type stars, and OB associations. We also estimate the probability 
of chance coincidences by using Monte Carlo simulations, and 
we provide a list of counterpart candidates of the gamma-ray 
sources. 



2. Cross-correlation of the First Fermi Catalog with 
massive young galactic objects 

There is not observational evidence of any YSO emitting 
gamma-rays so far. The new generation of gamma-ray tele- 
scopes, like Fermi, will make it possible to reach the necessary 
sensitivity level to detect these faint gamma-ray sources soon, in 
case the predictions were right. To identify those young objects 
that might be emitting gamma rays, we first took the recently 
published First Fermi Catalog ( |Abdo et al.||2010| l by the Fermi 
Collaboration and we excluded all known firm identifications, 
getting a list of 1392 sources. Then we crossed this list with cat- 
alogs of confirmed and well characterized YSOs and other type 
of young stars. We also did a Monte Carlo study to determine 
the probability of pure chance coincidences between the crossed 
catalogs. 

2.1. Catalogs 

The catalogs used in this study are listed in Table [T] Here we 
describe each one in more detail. 

The Fermi Large Area Telescope First Catalog ( |Abdo et"aL 
2010] ) contains the detected sources during the first 1 1 months of 
the science phase of the mission, which began on 2008 August 
4. This catalog contains 1451 gamma-ray sources detected 
and characterized in the 100 MeV to 100 GeV range with a 
typical position uncertainty of ~ 6'. After excluding the firm 
identifications from the original sample, we get 1392 sources. 
Most of them are located on the Galactic plane (see Figure[T]l. 



an arcsecond. 



The Red MSX Source (RMS) survey is an ongoing 
multi wavelength (from radio to infrared) observational program 
with the objective of providing a well-selected sample of 
MYSOs in the entire Galaxy (Urqu hart et al.| |2008] >. About 
~2000 MYSO candidates have been identified by comparing 
the colors of MSX and 2MASS point sources (at 8, 12, 14, and 
23 /im) with those of well known MYSOs. The survey also uses 
high-resolution radio continuum observations at 6 cm obtained 
with the VLA in the northern hemisphere and at 3.6 cm and 
6cm with ATCA in the southern hemisphere. They help to 
distinguish between genuine MYSOs and other types of objects, 
such as ultracompact HII regions, evolved stars, or planetary 
nebulae, which contaminate the sample. In addition to these 
targeted observations, archival data of a previous VLA survey of 
the inner Galaxy were used. This ongoing program has provided 
a sample of 637 well-identified MYSOs until now, which were 
used in our work. 



The Vllth catalog of Population I WR stars ( van der Hucht 
2001) contains 227 stars, with spectral types and bv photometry. 
In recent years, the number of WR stars has increased in 7 1 new 
stars, respect to the Vlth catalog and the coordinates have also 
been improved. The position uncertainty is close to a fraction of 



et al 



The catalog of Of-type stars is the one of Cruz-Gonzalez 
1974) 1, which contains 664 stars. The catalog provides 



m v , B-V, spectral type, radial velocity, radial component of the 
peculiar velocity, possible multiplicity of the object, and other 
characteristics for each source. The typical uncertainty in the 
star position is ~ 1'. 



Finally, the catalog of OB associations is the one by Mel' Nik 
& Efremov] ( |1995| >. This catalog contains 88 associations and 
provides distances to the association, number of stars, and size 
of the association along the Galactic latitude and longitude axes. 
The typical value of the size is ~ 20 - 30 pc. 

2.2. Spatial coincidences 

We crossed the Fermi catalog with the catalogs of young galactic 
objects mentioned above. We calculated the distance between 
two sources using the statistical parameter S ( Allington-Smith 
|et al.|1982| ): 



S = 



(Aacosd) 2 A6 2 



erf + <r 2 . 



where Aa and AS are the difference between the right ascension 
and the declination of the two compared sources, respectively, 
<r ab is the uncertainty in the position of the source, and (i, j) rep- 
resent the two sources. The error in the position of the Fermi 
sources is taken as the 95% confidence ellipse. The error in the 
position of the other compared sources is the precision in the 
coordinates for YSOs, WR stars and Of-type stars, and the an- 
gular size of the association in the case of OB associations. If 
S is lower than or equal to the unit, it means that the source 
position (YSO, WR, Of-type, or OB associations) is inside the 
95% uncertainty ellipse of the Fermi source, within its own posi- 
tion uncertainty, and that case is considered to be a coincidence. 
Massive YSO and protostar are point-like objects compared with 
the confidence contours of Fermi sources. The same is valid for 
WR and Of stars, either in binary systems or isolated. OB as- 
sociations are large systems that can contain more than a sin- 
gle gamma-ray source. We note that our study is based on two- 
dimensional coincidences, since we are comparing the equatorial 
coordinates of the sources. 



2.3. Monte Carlo analysis 

To determine the chance coincidences we used the Monte Carlo 
method for simulating sets of synthetic gamma-ray sources start- 
ing from the Fermi Large Area Telescope First Catalog. We fol- 
lowed a similar criteria to that used by Rome rcTet al.| ([1999) 
to search for the possible association of unidentified EGRET 
sources with other type of celestial objects. In this algorithm, 
the galactic coordinates of a gamma-ray source (/, b) are moved 
to new ones (l',b'). The new galactic longitude coordinate is cal- 
culated by doing I' = I +R\ x 360°, where Ri is a random number 
between and 1 that never repeats from source to source or from 
set to set. Since the distribution of Fermi sources is almost con- 
stant in Galactic longitude, we do not impose any constraint on 
this coordinate in the simulations. The sources in the Fermi cat- 
alog have a given distribution in galactic latitude (see Figure [TJ. 
In order to constrain the simulations with this distribution, the 
galactic latitude coordinate is calculated by doing V - b+RiV. 1°, 
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Table 2. Statistical results obtained from simulations. Latitude 
galactic coordinate has been constrained, while galactic longi- 
tude remains free. 



Galactic latitude (°) 



Fig. 1. Distribution of First Fermi Catalog sources in galactic 
latitude. 



where again R2 is a random number between and 1. Here, if the 
integer part of b' is greater than the integer part of b or the sign 
of b' is different than the sign of b, then V is replaced by b' — 1°. 



Table 1. List of the catalogs used in the study. 



Object 


Catalog 


#of 


type 




sources 


y-ray sources 


First Fermi Catalog" 


1392' 


YSO 


RMS Survey* 


637 


WR 


VII" 1 Catalog of Galactic Wolf-Rayet stars" 


227 


Of-type 


A catalog of galactic O stars and 


664 




the ionization of the low density 






interstellar medium by runaway stars d 




OB associations 


A new list of OB 


88 




associations in our Galaxy'' 





Notes. {a) 
1 2001k (</J 



Abdo et al. ( 2010 1, {b) Urquhart et al. ( 2008 ) 



if) 



van der Hucht 

Cruz-Gonzalez et al. \ \ 1 974) , w |Mej'Nik & Efremov| \\99$ 



{J> number of gamma-ray sources after excluding the firm identifications 



We simulated 1500 sets of synthetic Fermi sources and each 
set was compared with a fixed set of different kinds of objects: 
YSOs, WR stars, Of-type, and OB associations. In each simula- 
tion, we calculated the distance between the two sources using 
the statistical parameter S. 

For each kind of compared objects we calculated the aver- 
age number of coincidences and its standard deviation after all 
Monte Carlo simulations. We calculated the chance coincidence 
probability using the actual number of coincidences for each 
type of object and assuming a Gaussian distribution in the simu- 
lations. This probability allows us to know the reliability of our 
study. We also repeated this process moving the Fermi sources 
in 2°-bins in galactic latitude, i.e. replacing the galactic latitude 
coordinate by V — b + 7? 2 x 2°. If the integer part of V is greater 
than the integer part of b or the sign of V is different than the 
sign of b, then V is replaced by V - 2°. This binning allows us 
to keep the initial distribution in galactic latitude as well. 



Object 


Coincident 


Simulated 


Probability 


Simulated 


Probability 


type 


y-ray sources 


1° - bin 


1° - bin 


2° - bin 


2° - bin 


YSO 


12 


4.4±2.0 


1.8X10" 4 


3.6±1.8 


5.6x10"" 


WR 


2 


1.3±1.1 


2.9x10"' 


1 .2*1.1 


2.9x10"' 


Of-type 


5 


2.9±1.7 


1.1x10-' 


2.9±1.7 


1.1x10"' 


OB assoc. 


107 


72.5±8.0 


4.2xl0" 6 


72.8±8.0 


5.5xl0" 6 



3. Results 

The results from our statistical study are shown in Table [2] In 
this table we list, from left to right, the object type, the num- 
ber of coincidences between the each compared catalog and the 
original Fermi catalog, the simulated average number of coinci- 
dences, and the chance coincidence probability for each binning 
in galactic latitude. We find 12 gamma-ray sources spatially co- 
incident with YSOs, 2 with WR stars, 5 with Of-type stars and 
107 with OB associations. 

From the Monte Carlo analysis, we see that there is a strong 
correlation between gamma-ray sources and YSOs: the cata- 
log cross-check returns 12 coincidences between gamma-ray 
sources and YSOs. The Monte Carlo simulations, for the case 
of displacing the Fermi sources in l°-bins, returns an average of 
coincidences of 4.4±2.2 sources, which is the number of chance 
coincidences. This means that 7.6 of the 12 coincident Fermi 
sources (~ 63% of the total coincidences with a ~ 4cr confidence 
level ) should be associated with a probability of chance coinci- 
dence of 1.8xl0~ 4 . Similarly, for displacing the Fermi sources in 
2°-bins, we obtain an average of coincidences of 3.6+1.8. That 
result indicates that 8.4 of the 12 coincident sources (~ 70% of 
the total coincidences with a ~ 5cr confidence level) should be 
associated with a probability of 5.6 x 10~ 6 of being chance coin- 
cidences. In a similar way, the association of Fermi sources with 
WR and Of-type stars is unclear since the number of actual coin- 
cidences and the results of the Monte Carlo simulations are too 
similar and thus the probability of chance coincidence is too high 
(0.29 and 0.11 for WR stars and Of-type stars, respectively). 
The results for WR and Of-type stars are different from those 
obtained by Romero et al. ( 1999) for EGRET sources. The prob- 
ability of chance coincidences has increased while the number 
of candidates decreased in the case of WR stars. In the case of 
Of-type stars the probability of chance association has increased 
as the number of coincidences increased as well. The probability 
of chance coincidence with OB associations is as low as ~ 10 -6 . 
In this case, however, the nature of the gamma-ray emission is 
not clear, as it is discussed in section [331 



3.1. Young stellar objects 

The association of gamma-ray sources and massive YSOs has 
been suggested in the last years after studying a reasonable sce- 
nario for the production of non thermal emission ( |Araudo et al.| 
2007 ). However, this is the first time that the study of YSOs as 
gamma-ray sources is carried out in a statistical way, taking ad- 
vantage of the Fermi catalog. 

The results of our cross-check, shown in Table [3] indicate 
that 12 gamma-ray sources are positionally coincident with 23 
YSOs. In this table we present, from left to right, the Fermi 
source name, its J2000 equatorial coordinates, its positional un- 
certainty, the spectral y-ray index, the energy flux (E>100 MeV), 
the YSO name, its J2000 equatorial coordinates, the angular dis- 
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tance between the two compared sources, the distance to the 
YSO, its IR luminosity, and the mass of the star forming region 
where it is embedded. 

In what follows we present a case by case discussion of the 
gamma-ray fields. 

- 1FGL J054 1.1 +3542. This source is coincident with 
three YSOs: G173.6328+02.8604, G173.6339+02.8218 and 
G173.6882+02.7222. Their luminosities are below 4.8xlO 4 L . 
The three objects belong to the G173. 6036+02. 6237 complex 
(in the RMS Survey notation), which is situated at a distance of 
1.6 kpc. We also find two Herbig Haro like objects, GGD 5 and 
GGD 6 (Gyulbudag hian et al.|1978 i, within the error ellipse of 
the Fermi source. Outside the Fermi error box the complex also 
harbors three more YSOs and an HII region. 

- 1FGL J0647.3+0031. The YSO G212.0641-00.7395 is the 
only coincidence with this source. Its kinematic distance is 6.4 
kpc and it has a luminosity of 2.5xlO 4 L . This source belongs 
to the G21 1.9800-00.9710 complex, together with another YSO 
that lies outside the error box of the gamma-ray source. 

- 1FGL J1256.9-6337. We find G303.5990-00.6524 within 
the error ellipse of this gamma-ray source. It is a YSO with a 
bolometric luminosity of 8.3 x 1O 3 L and located at a kinematic 
distance of 11.3 kpc. It belongs to the G303.5670-00.6253 
complex that also harbors an HII region. 

- 1FGL J 13 15. 0-6235. This source is coincident with 
G305.4840+00.2248, which is a YSO with a luminosity of 
3.8 x 1O 3 L and located at a distance of 3.6 kpc. 

- 1FGL J1651.5-4602. The source G339.8838-01.2588 is 
coincident with this gamma-ray source. It is a YSO with a 
bolometric luminosity of 2.1xlO 4 L . It is located at 2.6 kpc 
from the Earth. This source has been detected in radio at 8.6 



GHz with an integrated flux of 2.6 mJy (Walsh et al. 



1998). 



There is no indication of whether the radio flux is non thermal 
or not. 

- 1FGL J1702.4-4147. This source is coincident with two 
YSOs: G344.4257+00451B and G344.4257+00451C, which 
form the G344.4120+00.0492 complex, together with two HII 
regions. They are located at 5 kpc from the Earth and both 
have a bolometric luminosity of 1.5xlO 4 L . These gamma-ray 

Dutra et al.|[2003l > 



source is also near to a cluster of stars (see 



and a molecular cloud (see Russeil & Castets 



( |2004D ). 



- 1FGL J1846.8-0233. We find a coincidence with the 
source G030.1981-00.1691. It has a bolometric luminosity of 
2.9xlO 4 L and is located at 7.4 kpc. Within the error ellipse of 
the Fermi source, there are several other sources, such as dark 
nebulae and HII regions. 

- 1FGL J1848. 1-0145. This source is coincident with 
two YSOs, G030.9726-00.1410 and G030.9959-00.0771, 
located at 5.7 kpc and with bolometric luminosities of 
3.9xlO 3 L and 5.1xlO 3 L , respectively. Both are part of the 
G03 1.145 1+00.0383 complex, which hosts five more YSOs 
(outside the Fermi error box), five diffuse HII regions, and ten 
HII regions. Within the error ellipse of the gamma-ray emission, 
there is an unidentified very high-energy gamma-ray source, 
HESS J 1848-01 8 ( [Chaves et aLl|2008} , which is probably the 
most suitable very high-energy candidate counterpart to the 



Fermi detection. 

- 1FGL J1853. 1+0032. This source has the biggest error 
ellipse (A#9g% ~ 0.5°). For that reason, there is a high number 
of sources within its location error box, including pulsars, 
supernova remnants and X-ray sources. The cross-match of 
the catalogs yields six coincidences. These six YSOs belong 
to three different complexes with different distance to the 
Earth: G032.8205-00.3300 and G033.3891+00.1989 belong 
to the G033. 1844-00.0572 complex, located at 5.1 kpc. The 
luminosities of these two sources are 1.7 and 1.1 xlO 3 L , 
respectively; G033.3933+00.0100 and G33.5237+00.0198 be- 
long to the G033. 6 106+00.0464 complex, located at a distance 
of 6.8 kpc and the luminosity is 7.9xlO 3 L for both sources; 
G034.0126-00.2832 and G034.0500-00.2977 belong to the 
G034. 03 13-00.2904 complex. The distance to these objects is 
13.3 kpc and their bolometric luminosities are 3.4xlO 4 L and 
2.4xlO 4 L , respectively. 

- 1FGL J1925.0+1720. Two YSOs are coincident with it: 
G052.2025+00.7217A and G052.2078+00.6890. The bolomet- 
ric luminosities are 1.5xlO 4 L and 2.0 xlO 4 L , respectively. 
They belong to the G052.2052+00.7053 complex which is 
located at 10.2 kpc, and hosts also a HII region. 

- 1FGL J1943. 4+2340. There is spatial coincidence with 
the source G059.783 1+00.0648, located at 2.2 kpc. It has 
a bolometric luminosity of 6.8xlO 4 L . This YSO has been 
detected in radio at 8.6 GHz with an integrated flux of 1.0 mJy 
( |Sridharan et al.|[2002| l. The infrared counterpart of this YSO 
is IRAS 19410+2336. It was observed by Chandra in 2002, 
which found hard X-ray emission from a number of sources 
within this high-mass star-forming region ( |Beuther eTa l. 2002a). 
The region has two cores where star formation takes place, 
with masses of 840 M and 190 M flSridharan et~aT[ [20021 
Beut her et al.|2002b) l. In the latter paper, it is proposed that the 
X-ray emission is produced by magnetic reconnection effects 
between the protostars and their accretion disks. The interaction 
of several molecular outflows, where the YSO from the RMS 
survey is located, and the combined effects of the stellar winds 
make a good scenario that might result in particle acceleration 
up to relativistic energies. 

- 1FGL J2040.0+4157. There is spatial coincidence with 
G08 1.5 168 +00. 1926. This YSO is located at 1.7 kpc and shows 
a bolometric luminosity of 7.O4xlO 2 L . There is a galaxy 
(2MASX J20395796+4159152) located at 2.3" from the posi- 
tion of that YSO. 



3.2. WR and Of-type stars 

We found two Fermi sources coincident with 15 WR stars in 
the Galactic center cluster and nine WR stars in the Quintuplet 
cluster. Some of these stars (see jvan der Huc ht 2006) show vari- 
ability and have been detected at X-rays with evidence of non- 
thermal emission. The results of our study are shown in Table|4] 
We cannot state that these two coincidences correspond to phys- 
ical associations given the high chance association probability 
obtained from the Monte Carlo study. The first one is in the di- 
rection of the galactic center, and there are several other sources 
that introduce confusion. The second one was suggested as po- 
tential association with the Pistol Star in |Abdo et al.| ( |20f0] l and 
it is situated in a very crowded field. 
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In the case of Of-type stars, which are the evolved state of 
O stars, and precursors of WR stars, our results show that five 
Fermi sources are coincident with five stars (see Table |5). The 
probability of chance coincidence is too high to state a phys- 
ical association. The sources 1FGL J 13 15. 0-6235 and 1FGL 
J 1853. 1+0032 also show positional coincidence with YSOs (see 
Table [3). The case of association with YSOs has a much lower 
value for the chance probability. The source 1FGL Jl 112.1- 
6041 is coincident with HD 97434, a multiple star system. 
1FGL J 13 15.0-6235 is located in a regions that harbors dark 
nebulae and molecular clouds. The Of star coincident with this 
source is HD 115071, which is a spectroscopic binary. Finally, 
1FGL J2004.7+3343 is coincident with HD 227465, and there 
is also the source G70.7+1.2 inside the error box of the Fermi 
detection, which might contain a Be star and an X-ray-emitting 
B star pulsar binary ( Cameron & Kulkarni|2007] l. 



3.3. OB associations 

Our study yields 107 Fermi sources positionally coincident with 
35 OB associations, which represent ~ 41% of the sample. 
We list the results in Table [6] available in the electronic ver- 
sion. We get this large number of gamma-ray sources because 
of the size of the OB associations. Most of them have angular 
sizes of ~ 1° in diameter or higher. Our results extend those of 
Rom ero et al.| (fl999 ), where they found 26 coincidences between 
EGRET sources and OB associations. Most of the OB associa- 
tions that they found are present in our results, along with other 
new candidates, as expected from the higher sensitivity of LAT. 
Although the number of sources has increased, the probability 
of chance associations is approximately the same as in |Romero| 
[etaLl ( fT999t . 

Most of the OB associations have fewer than five Fermi 
sources within their error boxes. There are seven OB associa- 
tions with five or more Fermi sources within their error boxes. 
There are four associations with a number of gamma-ray coin- 
cidences between five and ten (Ori 1 B, Ori 1 C, Car 2, and Cyg 
1,8,9), located at short distances from the Earth (less than 1 kpc) 
and located at galactic latitudes of ~ |15°| (association centroid 
position). There is an exception, Car 2, which is located at 2.2 
kpc and has a galactic latitude of -0.13°. The OB associations 
with the large number of gamma-ray coincidences (> 10) are 
those from Scorpius (Sco 2 A, Sco 2 B and Sco 2 D). Those as- 
sociations are located at ~ 170 pc on average and have very large 
angular sizes (~ 9.5° on average). 

There are five Fermi sources that are coincident with OB 
associations and YSOs at the same time: 1FGL J1256. 9-6337, 
1FGL J1315.5-6235, 1FGL J 1702.4-447, 1FGL J1943.4+2340, 
and 1FGL J2040.0+4157. In all cases the OB association over- 
laps both the Fermi source and the YSO In all cases but one, 
however, the distances to the YSO and the OB association are 
much too different. 

4. Discussion 

The fact that there are five Fermi sources coincident with both 
YSOs and OB associations makes us consider which is the 
chance probability of having coincidence of Fermi sources with 
YSOs alone. Using the Monte Carlo algorithm again, but taking 
this constrain into account, we obtain a mean value of coinci- 
dences of 2.8 + 1.7, and the chance probability for the seven 
Fermi sources coincident only with YSOs is ~ 1.6%. 

The brightest IR YSOs have more molecular mass avail- 
able for proton-proton collisions and Bremsstrahlung interac- 




34.5 35.0 35.5 

log(L-,) (ergs- 1 ) 

Fig. 2. IR luminosity versus gamma-ray luminosity above 100 
MeV of the brightest YSOs coincident with each Fermi source. 
The shaded area represents the 1-cr confidence interval of a least 
square fit taking eight data points randomly and repeating the fit 
one thousand times. It can be roughly seen that the higher the 
gamma-ray luminosity, the higher the IR luminosity. 



tions than those that are faint. It is expected then that the bright- 
est IR YSOs would show the highest gamma-ray luminosity. To 
test such a trend, we have plotted in Figure [2] the IR luminosity 
versus the gamma-ray luminosity of the brightest YSOs coin- 
cident with each Fermi source. The gamma-ray luminosity was 
calculated from the gamma-ray flux (see Table [3]) assuming that 
the distance to the gamma-ray source is equal to the distance 
to the corresponding YSO. We see that, under this assumption, 
there is a trend toward increasing the gamma-ray luminosity as 
the IR luminosity increases, too. This trend, however, is obtained 
using the whole data set, and according to the results from our 
Monte Carlo simulations, we should have four chance coinci- 
dences out of the 12 candidates. To see how this could affect 
the trend, we selected eight data points randomly and fitted them 
to a straight line by least squares. We repeated this process one 
thousand times, getting an average least square fit within one 
standard deviation. The limits of this fit are plotted in Figure 2. 
The increasing trend in the data is visible, although with a sig- 
nificant dispersion. This is not surprising considering the broad 
approach. 

The association of WR stars with gamma-ray sources has 
been discussed in the past by Kaul & Mitra ( 1997) and Romero 
|et al.| ( |T999] >. In both cases, the authors studied the positional co- 
incidence between WR stars and unidentified EGRET sources. 
In the later work, they find that two WR stars are of special in- 
terest: WR 140 and WR 142. The first one is a binary system, 
composed of a WC 7 plus an 04-5 star where the region of colli- 
sion of the winds seems to be a good place for particle accelera- 
tion and high energy emission. It should be mentioned, however, 
that WR 140 is a long-period binary with variability expected on 
time scales of years (Willia ms et al.|[T987] >. In the second one, 
WR 142, the hard X-ray emission and fast wind may indicate 
a colliding wind shock that could be explained by a companion 
close to this star ( Sokal et al. 2010| l. No companion, however, has 
been reported so far. In our work none of these mentioned WR 
stars appear to be coincident with any Fermi source. The poor 
statistical correlations found from the simulations does not allow 
us to be confident with any of the coincidences found. The case 
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of Of-type stars is also unclear since the probability of chance 
association is high. 

Finally, our results for the OB associations are not conclu- 
sive. We list the results in Table [6] available in the electronic 
version. The probability of chance coincidence is negligible 
(~ 10~ 6 ), but we do get several gamma-ray sources for each OB 
association. Thus, is very difficult to assign a specific counterpart 
to the gamma-ray emission. 

5. Conclusions 

We studied the two-dimensional coincidence between uniden- 
tified Fermi sources and catalogs of galactic young objects, 
such as YSOs, WR stars, Of stars, and OB associations. We 
have found a statistical correlation between gamma-ray sources 
and YSOs. The correlation with the early type stars with 
strong winds remain unclear, since the candidates are located 
in crowded fields with many other alternatives to the gamma-ray 
emission, and the probability of chance associations is high. In 
the case of OB associations, the probability of chance associa- 
tions is negligible. However, we cannot assign a specific coun- 
terpart to the gamma-ray emission because of the high angular 
size of most of OB associations. 

What we have presented here is the first statistical evidence 
for gamma-ray emission from massive YSOs. 
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Table 3. Positional coincidence of Fermi sources with MYSOs. 



Fermi Name 


RA 


Dec 


95% Semi 


Spectral Index r 


Flux(E>100MeV) 


MSX Name 


RA 


Dec 


Afl 


Distance* 




Mass 


(1FGL) 


o 


O 


Major Axis (°) 


(F oc E-' ) 


xl0~" erg cm" 2 s~' 




O 


o 


O 


(kpc) 


(xlOX) 


(Mo) 


J U J4 1 . 1 tjjtZ 


85 ^805 


35 7091 


1397 


2 41 +0 13 


1 6+0 5 


H173 ii37X+n"> Xllftd 




TJJ.OZUJJ 


12 


1 6" 


4 8 !/ 
















G173.6339+02.8218 


85.29592 


+35.83380 


(1.13 


1.6" 


3.2" 
















G173.6882+02.7222 


85.22758 


+35.73558 


0.05 


1.6" 






J0647.3+0031 


101.8417 


0.5289 


0.2150 


2.41±0.11 


1.9+0.5 


G212.0641-00.7395 


101.80567 


+0.43514 


0.10 


6.4* 


25' 




J1256.9-6337 


194.2474 


-63.6212 


0.1955 


2.26±0.12 


4.9+1.1 


G303.5990-00.6524 


194.35546 


-63.51650 


0.12 


11.3* 


8.3' 




j 1 3 1 5 0-6235 


198.7635 


-62.597 1 


0. 1 860 


2.3 1±0. 1 2 


6.9+0.0 


G305 .4840+00. 2248 


198.40016 


-62.53708 


] $ 


3.6* 


3.8/ 




J1651. 5-4602 


252.8831 


-46.0340 


0.2258 


2.21+0.07 


13.9+3.4 


G339. 8838-01 .2588' 


253.01942 


-46.14267 


0.14 


2.6'' 


21.0' 




J1702.4-4147 


255.6039 


-41.7859 


0.0800 


2.39+0.07 


8.7+2.0 


G344.4257+00.045 1 B 


255.53674 


-41.78303 


0.05 


5.0* 


15.0^ 
















G344.4257+00.0451C 


255.53587 


-41.78617 


0.05 


5.0'' 


15.0/ 




J1846. 8-0233 


281.7001 


-2.5628 


0.1262 


2.21+0.06 


9.3+2.3 


G030.1981-00.1691 


281.76274 


-2.51003 


0.08 


7.4* 


29.0' 




J1848. 1-0145 


282.0470 


-1.7605 


0.0859 


2.23+0.04 


9.5+3.2 


G030.9726-00.1410 


282.09178 


-1.80842 


0.07 


5.7' 


3.9^ 


1.9xl0 3 « 














G030.9959-00.0771 


282.04516 


-1,75808 


0.0044 


5.7* 


5.1' 


1.9X10 3 * 


J1853. 1+0032 


283.2884 


0.5366 


0.5207 


2.18+0.07 


5.7+1.7 


G032. 8205-00. 3300 


282.04436 


-1.75703 


0.34 


5.1* 


n.o> 
















G033.3891+00.1989 


282.89092 


+0.49750 


0.40 


5.1* 


11.0/ 
















GO33.3933+OO.O100 


283.06109 


+0.41528 


0.26 


6.8* 


7.9' 
















G033.5237+00.0198 


283.11179 


+0.53569 


0.34 


6.8* 


7.9' 
















G034.0126-00.2832 


283.60437 


+0.83239 


0.43 


13.3* 


34.0' 
















G034.0500-00.2977 


283.63454 


+0.85914 


0.47 


13.3* 


24.0' 




J1925.0+1720 


291.2748 


17.3485 


0.1443 


2.28+0.12 


2.4+1.0 


G052.2025+00.7217A 


291.24933 


+ 17.42169 


0.08 


10.2* 


15.V 
















G052.2078+00.6890 


291.28553 


+ 17.41317 


0.07 


10.2* 


20.0' 




J1943.4+2340 


295.8667 


23.6815 


0.1118 


2.23+0.11 


2.6+0.7 


G059.783 1+00.0648" 


295.79680 


+23.73433 


0.08 


2.2" 


6.8' 


840 and 190 * 


J2040.0+4157 


310.0154 


41.9533 


0.1970 


2.66+0.06 


7.9+1.2 


G081.5168+00.1926 


309.99066 


+41.98739 


0.04 


1.7' 


0.704' 





Notes. (1) Detected in radio at 8.6 GHz with integrated flux of 2.6 mjy; (2) Detected in radio at 8.6 GHz with integrated flux of 1.0 mjy. 
(3) Observed in the X-ray with Chandra (Beuther et al. 2002a); Distances: ia) distance to the complex, taken from the literature, <6) kinematic 
distance determined from the systemic velocity of the complex, (c) distance has been taken from the literature. Luminosities: {d) IRAS fluxes, 
{e} MSX 21 yum band flux using a scaling relationship determined from a comparison with sources where spectral energy distribution (SED) fits 
have been possible. W) SED fit to the available infrared fluxes (2MASS, MSX, MIPSGAL/1GA) and literature (sub) millimetre fluxes. w data obtained 
fr omhttp : //www.ast.le eds.ac.uk/cgi - bin/RMS/RMS_DATABASE.cgi. Masses: <s) from Rag an et al.||2006| l, m mass of the 2 cores identified 
in |Beuther etaLl ( |2002b> 



Table 4. Positional coincidence of Fermi sources with WR stars. 



Fermi Name 


RA 


Dec 


95% Semi 


r 


Flux(E>100MeV) 


Star 


RA 


Dec 


A0 


Distance 


(1FGL) 


(°) 


(°) 


Major Axis (°) 


(F oc E~ Y ) 


Xl0" 10 erg cm" 2 s" 1 




o 


O 


(°) 


(kpc) 


J 1745. 6-2900 


266.420 


-29.014 


0.019 


2.26±0.03 


7.1±0.7 


WR 101a 1 


266.41454 


-29.00950 


0.01 


8.0 


J 1746.4-2849 


266.618 


-28.818 


0.10 


2.2±0.5 


4.6±0.0 


WR 102c 2 


266.54667 


-28.81822 


0.06 


8.0 



Notes. (1) WR stars 101 b, c, d, e, f, g, h, i, j, k, 1, m, n, o, belonging to the Galactic Center Cluster lie at distances less than 20" and are also within 
the error box of the gamma-ray source; <2) WR stars 102 d, e, f, g, h, i, j, k belonging to the Quintuplet Cluster lie at distances less than 1' and are 
also coincident with the gamma-ray source. 



Table 5. Positional coincidence of Fermi sources with Of stars. 



Fermi Name 


RA 


Dec 


95% Semi 


r 


Flux(E>100 MeV) 


Star 


RA 


Dec 


A6» 


Distance 1 


(1FGL) 


o 


C) 


Major Axis (°) 


(F oc E~ r ) 


xl0~ 12 erg cm" 2 s" 1 




n 


(°) 


O 


(kpc) 


J0005. 1+6829 


1.2841 


68.4883 


0.4427 


2.58±0.12 


1.7±0.5 


BD+67 1598 


1.225 


68.167 


0.32 


1.07 


Jl 112.1-6041 


168.0486 


-60.6929 


0.0461 


2.12±0.05 


14.1 + 1.6 


HD 97434 


167.975 


-60.683 


0.04 


2.67 


J1315.0-6235 


198.7644 


-62.5971 


0.1860 


2.31±0.12 


6.9±0.0 


HD 115071 


199.000 


-62.583 


0.11 


0.74 


J1853.1+0032 


283.2887 


0.5369 


0.5207 


2.18±0.07 


5.7±1.7 


BD-0 3584 


283.400 


0.567 


0.12 


2.18 


J2004.7+3343 


301.1855 


33.7171 


0.1320 


2.28±0.08 


5.2±0.8 


HD 227465 


301.125 


33.700 


0.05 


3.48 



Notes. (1) distances from |Cruz-Gonzalez et al.|{T974) 
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